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Steven  A.  Orszag,  Principal  Investigator 
Department  of  Mechanical  and  Aerospace  Engineering 
Princeton  University 
Princeton,  NJ  08544 

In  the  attached  papers,  we  summarize  work  done  on  this  research 

project.  The  major  results  include: 

1.  Development  and  application  of  the  renormalization  group 
method  to  the  calculation  of  fundamental  constants  of  turbulence, 
the  construction  of  turbulence  transport  models,  and  large-eddy 
simulations. 

2.  The  application  of  RNG  methods  to  turbulent  heat  transfer 
through  the  entire  range  of  experimentally  accessible  Reynolds 
numbers. 

3.  The  discovery  that  high  Reynolds  number  turbulent  flows  tend 
to  act  as  if  they  had  weak  nonlinearities,  at  least  when  viewed  in 
terms  of  suitable  ’quasi-particles.’  This  discovery  suggests  that  tur¬ 
bulence  acts  as  if  there  were  a  significant  scale  separation  in  the 
flows,  even  though  turbulence  does  not  appear  to  have  such  scale 
separation.  These  ideas  seem  to  provide  some  justification  for 
eddy  transport  ideas  that  have  proven  so  useful  in  engineering 
descriptions  of  turbulence. 

4.  The  further  analysis  of  secondary  instability  mechanisms  in 
free  shear  flows,  including  the  role  of  these  instabilities  in  chaotic, 
three-dimensional  free  shear  flows. 


The  further  development  of  numerical  simulations  of  turbulent 
spots  in  wall  bounded  shear  flows. 

The  study  of  cellular  automata  for  the  solution  of  fluid 
mechanical  problems.  In  particular,  a  realistic  assessment  of  the 
utility  of  these  new  methods  for  complex  high  Reynolds  number 
flow  problems  has  been  given. 

The  clarification  of  the  relationship  between  the  hyperscale  ins¬ 
tability  of  anisotropic  small-scale  flow  structures  to  long- 
wavelength  perturbations  and  the  cellular  automaton  description 
of  fluids. 

The  development  of  efficient  methods  to  analyze  the  structure 
of  strange  attractors  in  the  description  of  dynamical  systems.  This 
includes  both  the  computation  of  Lyapunov  exponents  and  the 
computation  of  dimensions  of  attractors. 

The  analysis  of  hyperscale  instability  as  a  mechanism  for  desta¬ 
bilization  of  coherent  flow  structures. 

Further  details  are  given  in  the  attached  papers. 
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Analogy  between  hyperscale  transport  and  cellular  automaton  fluid 
dynamics 

Victor  Yakhot,  Bruce  J.  Bayly,  and  Steven  A.  Orszag 

Applied  and  Computational  Mathematics,  Princeton  l/niversity,  Princeton,  New  Jersey  08544 
(Received  24  February  1986;  accepted  16  April  1986) 

It  is  argued  that  the  dynamics  of  a  very  large  scale  ( hyperscale )  flow  superposed  on  the  stationary 
small-scale  flow  maintained  by  a  force  f(x)  is  analogous  to  the  cellular  automaton 
hydrodynamics  on  a  lattice  having  the  same  spatial  symmetry  as  the  force  f. 


While  real  fluids  consist  of  discrete  particles,  they  can  be 
regarded  as  continuous  media  at  scales  that  are  much  larger 
than  the  typical  intermolecular  distance,  and,  on  these 
scales,  they  can  be  described  by  the  equations  of  continuum 
hydrodynamics.  These  equations  are  quite  insensitive  to  the 
details  of  the  molecular  dynamics;  the  microscopic  interac¬ 
tions  affect  only  the  viscosity  coefficient.  Microscopically 
dissimilar  fluids  can  be  described  by  the  Navier-Stokes 
equations,  although  the  microscopic  properties  of  different 
fluids  may  be  reflected  in  a  very  wide  range  of  viscosity  coef¬ 
ficients. 

The  lack  of  dependence  of  the  hydrodynamics  on  the 
microscopic  properties  of  the  fluids  is  the  basis  for  the  recent 
interest  in  discrete  approximations  to  molecular  dynamics 
or  Cellular  automata  (CA’s).'“*  Cellular  automata  are  dis¬ 
cretely  and  locally  linked,  finite  state  machines.  The  "mole¬ 
cules”  in  a  CA  fluid  move  in  discrete  steps  over  the  lattice 
sites  and  interact  according  to  a  well  defined  set  of  rules  that 
typically  conserve  momentum  and  the  total  number  of  parti¬ 
cles.  The  hydrodynamic  behavior  of  the  CA  fluid  is  given  by 
the  evolution  of  the  average  macroscopic  properties  of  the 
system  (“slow"  modes).  Some  limitations  of  CA  hydrodyn¬ 
amics  have  been  discussed  in  Ref.  5. 

The  first  lattice  model  of  a  fluid  was  introduced  by  Har¬ 
dy,  de  Passis,  and  Pomeau  (HPP).'  Recently,  new  models, 
which  are  modifications  of  the  HPP  ideas,  have  led  to  simu¬ 
lations  of  two-dimensional  fluid  motions  that  appear  to  be 
compatible  with  experimental  data.2-*  Hydrodynamic  equa¬ 
tions  for  such  a  CA  fluid  can  be  derived  using  techniques 
based  on  the  Chapman-Enskog  expansion,  as  in  the  kinetic 
theory  of  gases.  It  has  been  found  that  the  form  of  the  contin¬ 
uum  equations  for  a  given  CA  fluid  depends  strongly  on  the 
symmetry  properties  of  the  lattice.  In  particular,  the  HPP 
lattice  gas  model  based  on  the  two-dimensional  square  lat¬ 
tice  leads  to  anisotropic  viscosity  and  anisotropic  nonlinear 
terms  in  the  resulting  continuum  dynamics.  However,  a  re¬ 
gular  hexagonal  lattice,  again  in  two  dimensions,  introduced 
by  Frisch.  Hasslacher.  and  Pomeau2  is  symmetric  enough  to 
produce  the  Navier-Stokes  equations  with  isotropic  viscos¬ 
ity  and  nonlinear  terms  It  has  been  pointed  out  by  Wolfram* 


that  in  three  dimensions  none  of  the  space-filling  crystallo¬ 
graphic  lattices  have  sufficient  symmetry  to  guarantee  the 
isotropy  of  the  corresonding  hydrodynamic  equations. 
However,  icosahedral  symmetry  would  produce  isotropic 
viscosity  and  nonlinear  terms.6  Unfortunately,  no  periodic 
lattice  has  such  symmetry  (with  the  possible  exception  of 
some  recently  conjectured  quasicrystal  structures).7 

Hydrodynamic  equations  are  derived  from  the  micro¬ 
scopic  equations  of  motion  by  averaging  over  small  scales.  It 
is  natural  to  pose  the  following  problem:  Let  us  consider  a 
viscous  fluid  driven  by  a  force  which  generates  a  stationary 
field  vf  on  the  small-scale  /.  The  equation  of  motion  for  the 
perturbation  ch  defined  on  scales  L  which  are  much  larger 
than  the  scale  l  of  the  basic  flow  can  be  derived  by  averaging 
over  the  small-scale  velocity  field  v‘ .  The  resulting  equations 
describe  hyperscale  hydrodynamics.  To  see  that  this  problem 
arises  quite  naturally,  let  us  imagine  a  system  of  microscopic 
particles  driven  by  an  external  force  f.  The  molecules  partici¬ 
pate  in  two  kinds  of  motion,  one  related  to  thermodynamic 
noise  and  the  other  caused  by  the  external  force.  Filtering 
out  the  smallest  (thermodynamic)  scales,  one  derives  the 
Navier-Stokes  equation  subject  to  the  external  force.  If  we 
also  average  over  the  scales  corresponding  to  the  force  f,  the 
resulting  equation  will  not  necessarily  be  the  Navier-Stokes 
equation,  but  will  rather  be  an  equation  describing  the  large- 
scale  (hyperscale)  motion  that  does  not  explicitly  include 
the  external  force. 

Some  examples  of  hyperscale  dynamics  have  been  con¬ 
sidered  in  Ref.  8.  It  has  been  shown  that  a  system  of  square 
vortices  (eddies)  gives  rise  to  the  equation  of  motion  for  the 
velocity  perturbation  at  large  scales  with  an  anisotropic  vis¬ 
cosity: 

v,  =  v(  1  -t-  j  Re2  -  1  Re2  sin  2d>),  ( 1 ) 

where  vt  is  the  effective  viscosity  in  the  direction  (cos<£, 
sin  <b).  In  this  and  all  subsequent  equations.  Re  denotes  the 
Reynolds  number  of  the  small-scale  flow;  the  formulas 
quoted  are  the  lowest-order  nontrivial  results  of  a  hierarchy 
of  successive  smoothing  approximations."  A  plane  parallel 
system  of  eddies  also  leads  to  an  anisotropic  viscosity  "  1  low- 
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ew.  a  system  of  triangular  vortices  having  hexagonal  sym¬ 
metry  (invariance  under  rotation  by  60°)  generates  an  iso¬ 
tropic  viscosity  coefficient  for  the  hyperscale  motion: 

v,  =  v(l  +!  Re2).  (2) 

The  analogy  with  cellular  automata  is  striking:  in  two 
dimensions  only  the  triangular  lattice  and  triangular  set  of 
vortices  produce  an  isotropic  equation  for  the  large-scale 
velocity  fluctuations.  Moreover,  it  has  been  shown  by  Siva- 
shinsky10  that  hyperscale  hydrodynamics  is  not  Galilean  in¬ 
variant  because  the  averaged  Navier-Stokes  equation  with 
the  forcing  term  is  not.  The  same  holds  for  the  CA  hydro¬ 
dynamics:  It  has  been  pointed  out2  that  the  continuum  equa¬ 
tions  following  from  cellular  automaton  models  are  not  Ga¬ 
lilean  invariant  as  a  result  of  the  discrete  lattice  underlying 
the  model. 

The  dynamics  of  hyperscale  flows  superposed  on  small- 
scale  flows  in  three  dimensions  have  been  studied  in  Refs.  9 
and  11.  The  analogs  of  steady  cellular  flows  in  two  dimen¬ 
sions  are  the  family  of  so-called  Beltrami  flows  in  three  di¬ 
mensions,  defined  by  finite  Fourier  sums  of  the  form 

v(x)  =  ^j\/4(Q)(n  +  fQXn)e'°'*/l.  (3) 

Here  S’  is  a  finite  set  of  unit-magnitude  vectors  Q,  n  is  a  unit 
vector  perpendicular  to  Q,  and  the  complex  amplitudes 
A(Q)  satisfy  the  reality  condition  A(  —  Q)  —A  *(Q). 
These  flows,  like  two-dimensional  cellular  flows,  are  exact 
steady  solutions  of  the  inviscid  fluid  equations,  and  can  be 
maintained  in  viscous  fluid  by  the  action  of  an  externally 
imposed  body  force. 

The  form  of  the  equation  of  motion  for  hyperscale  flow 
on  the  small-scale  flow  (3)  involves  the  fourth-rank  ten¬ 
sor9-'* 


AW  =  J>(Q)|J(^„  ~  2 StQjQ, 

—  <5  jiQ,Qk  —  26  tJ  Qk  Qi  +4  Q,Q  jQkQi), 
which  is  isotropic  only  if 

T\A(Q)\iQlQj  =A6fi, 
o 

2  =^(<5y(5k,  •+•  SitSkj  +  6lkSji), 

OKS 

for  some  constants  A  and  fi.  Clearly,  the  more  vectors  Q  we 
have  in  S,  the  better  our  chances  of  being  able  to  select  ampli¬ 
tudes  A(  Q )  of  the  corresponding  components  so  as  to  obtain 
an  isotropic  tensor  N.  We  shall  give  some  examples  to  illus¬ 
trate  the  connection  between  the  underlying  small-scale  flow 
structures  and  the  resulting  hyperscale  dynamics. 

The  simplest  Beltrami  flow  in  the  family  (3)  has  only 
two  Fourier  components  with  the  wave  victors 
±  Qo  =  (  ±  1.0,0)  and  amplitudes  A(  ±  Q„)  =  c/2;  v' 
=  c  (O.cos  x/l,  —  sin  x/I),  A  hyperscale  perturbation  with 
the  wave  vector  (O.k.O)  then  obeys  the  equation  of  motion 

f,  =  —  vA:  :( 1  —  }  Re:)u„  i  .=;0. 
f,  =  -  vk  2(  1  +  j  ReJ)i:v 

The  effective  viscosity  for  the  r,  component  is  negative,  and 


the  hyperscale  flow  is  therefore  unstable,  if  the  small-scale 
Reynolds  number  Re=c/  / v  exceeds  ^2. 

The  so-called  ABC  flow  is  obtained  when  S  consists  of 
six  wave  vectors  located  at  the  vertices  of  an  octahedron: 


5=  {Q}  = 


(4) 


with  all  amplitudes  having  the  same  modulus  i^(Q)|  =  c/ 
\T2,  where  c  is  again  the  rms  fluid  velocity.  Now  the  hyper¬ 
scale  equation  takes  the  form 


4~va=  -  vAt2[<5a6  +  Re1  J/„*(k)]iv, 
at 

where  M  is  a  matrix  that  depends  only  on  the  unit  vector  k  in 
the  direction  of  k.  The  eigenvalues  of  M  are  all  greater  than 
or  equal  to  zero,  with  equality  occurring  only  if  k  lies  in  one 
of  the  coordinate  planes.  The  hyperscale  flow  is  therefore 
stable  although  still  somewhat  anisotropic. 

A  more  complex  flow  with  similar  structure  to  the  ABC 
flow  can  be  obtained  by  augmenting  the  wave  vector  set  by 
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and  assigning  the  corresponding  Fourier  modes  amplitudes 
A(Q)  =A(c/yJT2),  QeS1.  Here  A  is  a  parameter:  setting 
A  **  0  recovers  the  ABC  flow  with  its  anisotropic  equation  of 
motion,  but  if  A  is  raised  to  the  special  value  of  9/(  156) 1/2 
then  the  conditions  for  isotropy  of  the  tensor  Nyk/  are  satis¬ 
fied,  and  we  obtain  an  isotropic  equation  of  motion  for  the 
hyperscale  modes.  This  example  illustrates  the  fact  that  a 
small-scale  flow  can  be  constructed  to  have  different  large- 
scale  properties  from  a  discrete  lattice  gas  with  the  same 
spatial  symmetry  group. 

Going  on  to  more  and  more  symmetric  flows,  it  turns 
out  that  icosahedral  or  dodecahedral  symmetry  in  the  small- 
scale  flow  gives  exact  isotropy  to  the  hyperscale  dynamics. 
For  example,  the  wave  vector  set  for  the  icosahedral  flow  is 


where  r  is  the  golden  ratio  (1  +  v5)/2  and  j  =  (5  +  \  5)/2. 
It  is  easily  checked  that  the  tensor  N„;,  is  isotropic  for  this 
flow,  provided  that  the  amplitudes  of  the  modes  are  all  cho¬ 
sen  equal.  The  hyperscale  properties  of  the  icosahedral  flow 
appear  to  be  indistinguishable  from  those  of  the  exactly  iso¬ 
tropic  flow  obtained  as  the  limiting  case  of  flows  with  more 
and  more  Fourier  components  distributed  uniformly  on  the 
unit  sphere. 

It  is  interesting  to  observe  that  not  only  docs  the  mo¬ 
mentum  equation  for  the  hyperscale  modes  take  the  classical 
form,  but  so  does  the  hyperscalc  diffusion  equation.  The  ten¬ 
sor  that  enters  the  correction  for  the  effective  diflusiuty  in 
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(he  lowest  smoothing  approximation’"'  is  (he  sum  of  the 
“projection  tensors" 

7*.,=  V^-G.o,). 

O'  s 

which  is  isotropic  for  all  the  aforementioned  flows  except  the 
simple  flow  with  only  two  Fourier  components.  Indeed,  it 
has  been  demonstrated  independently  of  the  smoothing  the¬ 
ory  that  a  passive  contaminant  disperses  diffusively  in  the 
icosahedral  and  augmented  cubic  flows.  Simulations  of  par- 
tide  dispersion  in  these  flows  demonstrate  that,  for  large 
times,  almost  all  particles  migrate  away  from  their  starting 
point  with  a  finite  effective  diffusivity. 

The  analogy  between  the  CA  and  hyperscale  hydrody¬ 
namic  descriptions  of  the  fluids  goes  even  further.  The  equa¬ 
tions  for  the  large-scale  velocity  field  derived  in  Refs.  7-10 
are  based  on  the  neglect  of  the  higher-order  nonlinear  terms 
generated  by  the  scale  elimination  procedure  and  thus  they 
are  valid  when  the  ratio  u1-  <vl.  The  same  holds  for  the  CA 
hydrodynamics:3  the  Navier-Stokes  equation  is  an  approxi¬ 
mation  valid  only  when  the  Mach  number  Ma  =  u/v,h  <1, 
where  v,h  is  the  velocity  of  the  particles  on  the  lattice. 

Based  on  the  analogy  between  hyperscale  hydrodyna¬ 
mics  and  the  CA  description  of  the  fluids,  we  argue  that 


lattice  gas  models  are  equivalent  to  the  Navier-Stokes  equa¬ 
tion  with  an  external  force  having  the  symmetry  of  the  lat¬ 
tice. 
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a  canal  with  a  porous  plate  located  at  a  step 
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A  linear  shallow  water  theory  of  waves  inlt^anal  is  considered/The  energy  dissipated  by  a  porous 
plate  located  at  a  step  where  both  the  depth  ahdwidth  of  the  Oanal  changes  abruptly  is  calculated. 
It  is  found  that  a  broadening  of  the  canal,  combined  with  /plate  with  optimal  properties,  is  an 
efficient  way  of  dissipating  wave  energy.  It  is  also  pHf  sibl^/j  have  a  strong  wave  damping  without 
any  reflection. 


Lamb1  gave  the  reflection  and  transmission  coefficients 
for  long  waves  at  a  step  in  a  canal.  Both  the  width  and  d^pth 
of  the  canal  may  change  abruptly  at  the  step.  The  detailed 
analysis  of  Bartholomeusz3  verified  that  Lamb’s  solution  is 
correct  to  the  leading  order  in  a  long-wave  expulsion  (see 
also  Mei3).  We  will  generalize  Lamb’s  solution  to  the  case  of 
a  thin  porous  plate  located  at  the  step,  and  dricuss  the  dissi¬ 
pation  of  wave  energy.  Thereby  the  recent  york  of  Chwang 
and  Dong4  is  extended. 

We  consider  a  straight  horizonurf  canal  of  uniform 
depth  A,  and  width  b,  for*  <0,  where**  is  a  coordinate  along 
the  canal.  For  x  >  0  the  depth  is  hA nd  the  width  b2.  At  the 
step  ( x  =  0)  there  is  a  thin  vertkal  porous  plate  of  thickness 
d  and  permeability  K.  The  kinematic  viscosity  of  the  fluid  is 
v,  but  the  flow  outside  the  plate  is  assumed  inviscid.  The  flow 
inside  the  porous  plate  is  assumed  to  be  governed  by  Darcy’s 
law,  which  may  be  stated  as 

u  =  (Kg/%-dHyt  -  y:) ,  (1) 


valid  in  the  long-wave  limit.  Here#  is  the  gravitational  accel¬ 
eration,  y,  and  are  the  surface  elevations  on  each  side  of 
the  plate,  and  u  is  the  average  horizontal  velocity  in  the  sec¬ 
tion  covered  by  the  plate: 

u  =\|/  min(  l,42/b|)min(l,/iJ//t,) 

=  M\min(/>1/6;,l)min(A,//i2,l)  .  (2) 

Here  u,  and  t\are  the  average  horizontal  velocities  next  to 
the  plate,  for  x  <,0  and  jr  >  0,  respectively. 

Following  Land),1  we  calculate  the  transmission  coeffi¬ 
cient  T,  and  the  refl«rion  coefficient  R ,  for  waves  incident 
from  x  <  0-. 

T,  =  2/(G  ~  '  -c  1  4- 
R,  =  (G-1  4  1  -  r)/(c\^l  4-  r)  . 

In  F.qs.  (3)  and  (4)  we  have  introduced  the  dinieiv 
parameters 

G  ~  (K  / vd)c .  tnin(/> ,//>.,! ) mm (6, /A., I ) 


(3) 

(4) 

'.unless 

(  5  ) 
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Reynolds  Number  Scaling  of  Cellular-Automaton  Hydrodynamics 
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We  argue  that  the  computational  requirements  for  presently  envisaged  cellular-automaton  simu¬ 
lations  of  continuum  fluid  dynamics  are  much  more  severe  than  for  solution  of  the  continuum 
equations. 

PACS  numbers:  47.10. +  g 


It  has  recently  been  suggested1,2  that  cellular  auto¬ 
mata  (CAs)  [defined  as  discretely  and  locally  linked, 
finite-  (and  few-)  state  machines]  may  be  an  effective 
way  to  compute  complex  lluid  flows.  These  automata 
have  the  advantage  that  they  may  be  simply  and 
perhaps  inexpensively  constructed  with  use  of  specially 
designed  parallel  hardware.  With  suitable  interaction 
rules,  it  has  been  argued,1, 2  the  space-time  average 
kinetic  behavior  of  the  CA  system  follows  the  in¬ 
compressible  Navier-Stokes  dynamical  equations. 
While  the  Navier-Stokes  equations  for  continuum 
fluids  can  be  calculated  efficiently  on  parallel- 
architecture  machines,  it  is  probably  easier  to  make  ef¬ 
ficient  use  of  the  parallel  architecture  with  CAs.  In 
this  Letter,  we  wish  to  point  out  that  there  are  some 
considerations  that  require  resolution  before  these 
methods  can  be  considered  to  be  a  viable  alternative  to 
traditional  continuum  mechanical  methods  for  high— 
Reynolds-number  fluid  dynamics. 

Let  us  compare  the  resolution  and  work  require¬ 
ments  for  a  CA  simulation  of  a  high- Reynolds- 
number  flow  with  those  of  direct  numerical  solution  of 
the  incompressible  Navier-Stokes  equations.  It  is  well 
known3,4  that,  at  Reynolds  number  jVRe,  the  Kolmo¬ 
gorov  and  Batchelor-Kraichnan  theories  of  three-  and 
two-dimensional  equilibrium  range  dynamics,  respec¬ 
tively,  predict  that  the  range  of  excited  scales  is  of  or¬ 
der  ,V^4  and  Nl'1  and  the  computational  work  re¬ 
quired  to  calculate  a  significant  time  in  the  evolution 
of  large-scale  flow  structures  is  of  order  .V^e  and  <V^eJ 
in  three  and  two  dimensions,  respectively. 

The  suggested  evolution  rules  for  CAs  to  reproduce 
hydrodynamic  behavior  are  based  on  conservation  laws 
of  mass,  momentum,  and  energy.  Dissipation  is 
modeled  through  the  thermalization  of  coherent  hy¬ 
drodynamic  modes.  Therefore,  the  lattice  resolution 
of  the  CA  calculation  must  be  much  finer  than  that  of 
the  hydrodynamic  simulation,  the  latter  requiring  the 
retention  of  only  those  degrees  of  freedom  describing 
motions  on  scales  of  the  dissipation  range  or  larger. 
Thus,  the  lattice  spacing  a  must  be  smaller  than  the 
dissipation  scale  rj  in  the  turbulent  fluid. 

We  now  discuss  some  conditions  that  CA  models 
should  satisfy  to  describe  high- Reynolds-number  fluid 


flows.  We  present  three  successively  more  restrictive 
arguments  that  show  that  rj/a  must  grow  rapidly  with 
Reynolds  number. 

Signal-to-notse  ratio. — The  hydrodynamic  velocity  in 
the  CA  simulation  is  calculated  by  subdividing  the 
computational  domain  into  cells  with  linear  dimen¬ 
sions  »  a,  averaging  over  the  CAs  within  a  (finite) 
cell,  and  smoothing  (filtering)  the  resulting  (noisy) 
velocity  field.  Thus,  the  hydrodynamic  velocity  at  a 
point  x  is  the  (space-lime)  filtered  velocity  of  the  CAs 
in  the  cell  Cx  centered  at  x,  vH(x)  -  {v(x)),  where 
the  local  velocity  in  Cx  is 


where  n  is  the  number  of  occupied  sites  /  within  the 
cell.  We  assume  that  the  possible  velocity  values  at  an 
occupied  CA  site  are  v,  -  ±  v,h  where  v,h  is  the  con¬ 
stant  (thermal)  velocity  over  the  CA  grid.  At  low 
Mach  numbers,  vH  «  t/,„.  In  this  case,  the  fluctua¬ 
tions  in  v(x)  are  of  order  rr~l/JvIh.  In  order  that  the 
hydrodynamic  velocity  found  in  this  way  may  be  a 
good  representation  of  the  continuum  hydrodynamics, 
it  is  necessary  that  the  noise  n_1/:v,h  be  small  com¬ 
pared  to  the  smallest  significant  hydrodynamic  veloci¬ 
ty.  The  smallest  significant  hydrodynamic  velocity  is 
the  eddy  velocity  on  scales  of  order  of  the  dissipation 
scale  tj.  In  three  dimensions.  ■n~0((t/v})~u*)  and 
the  eddy  velocity  on  the  scale  of  77  is  vv  -  0  ( ( « i- ) 1/4 ) . 
Here  v  is  the  viscosity  and  c  is  the  turbulent  energy 
dissipation  rate  per  unit  mass.  Thus,  we  require  that 
the  number  of  CA  sites  n  within  a  cell  of  size  17  be  at 
least 

n»vjltv)m.  (2) 

Since  «-0((/3/L)  where  U  is  the  large-scale  rms 
fluctuating  velocity  and  L  is  the  associated  large-scale 
length  of  these  velocity  fluctuations,  we  find  that 
n  »  where  ,VRe  —  UL/v  is  the  Reynolds 

number  and  V/  -  U! i>lh  is  the  Mach  number.'  Since 
the  number  of  cells  of  size  rj  wuhm  a  three- 
dimensional  turbulent  eddy  of  size  I.  scales  as  .VRj\ 
the  overall  number  of  CA  sites  must  increase  at  least 
as  .VJi'Wlf1. 
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Since  the  effective  evolution  time  of  the  fluid  sys¬ 
tem  is  L/  L\  while  the  lime  step  on  the  CA  lattice  is 
a/vtn,  it  follows  that  the  CA  simulation  icquires  at 
least  L/a\l  steps  in  time.  Since  the  computational 
work  for  each  site  update  is  of  order  1,  it  follows  that 
the  CA  simulation  requires  at  least  of  order  (.VRi./ 
At  )ll/J  work. 

In  summary,  the  above  signai-to-noise  considera¬ 
tions  suggest  the  following  lower-bound  estimates  for 
the  computer  storage  S  and  work  W  for  CA  simula¬ 
tions  of  high-Reynolds  number,  low-Mach-number 
flows  (where,  for  reference,  we  include  the  corre¬ 
sponding  estimates  for  the  continuum  Navier-Siokes 
equations):  forCA  (2D), 

S-N&/M2. 

for  Navier-Siokes  (2D),6 

s-/vRe.  M'-Atf,*; 

forCA  (3D), 

S-zV^/A/2,  W  -  (iVRe/A/  )11/3; 
for  Navier-Stokes  (3D),7 

S-.V6/e4.  W-Nlr 

Upper  bound  for  the  Reynolds  number. — A  more 
stringent  condition  on  the  Reynolds-number  depen¬ 
dence  of  the  minimal  number  of  lattice  sites  in  a  CA 
simulation  of  hydrodynamics  is  found  as  follows.  If 
the  discrete  velocity  of  the  CAs  is  z  (again,  the 
thermal  velocity  or  sound  speed  on  the  CA  lattice)  and 
the  lattice  spacing  is  a ,  then  the  kinematic  viscosity  v 
on  the  lattice  is  at  least  of  order  va.  For  the  CA  to 
give  a  self-consistent  hydrodynamic  simulation,  the 
viscosity  determined  on  the  “molecular”  level  must 
equal  the  viscosity  governing  the  dissipation  of  the  hy¬ 
drodynamic  modes.  Thus  the  Reynolds  number  of  the 
simulated  fluid  can  be  at  most  L'L/v  or  ML/a.  Since 
the  number  .V  of  CA  sites  in  the  lattice  is  of  order 
(L/a  )d ,  where  J  is  the  dimension  of  space,  we  obtain 
the  result  that  .V  must  be  at  least  of  order  (.VRe,  ,V/  )J. 
As  above,  the  CA  simulation  of  the  flow  requires  at 
least  L/aAl  steps  in  time.  It  follows  that  the  CA  simu¬ 
lation  requires  at  least  of  order  i.\’KJ.\f)d  memory 
and  of  order  ,VR*  1  /.\tJ*2  work. 

These  estimates  for  storage  5  and  work  IF  based  on 
lower  bounds  for  the  effective  viscosity  on  the  lattice 
are  of  order 

S-<.VKt/.U):.  W-Ni'/At* 
forCA  CD). 

S-  A'u,..  W  -  VRe- 
for  Navier-Siokes  CD), 

s-  <  vKe.  ,vn\  h  -yjt,A/' 
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lor  C  A  (31)). 

S  =  .V^.  M'-.Vi, 

for  Navier-Siokes  (3D). 

Hydrodynamic  Jluauatiuns. — The  CA  system  will 
yield  a  self-consistent  continuum  hydrodynamic  de¬ 
scription  only  if  the  thermal  energy  fluctuations  on  hy¬ 
drodynamic  spatial  scales  are  small  compared  to  the 
energy  of  the  hydrodynamic  modes  on  corresponding 
length  scales.  If  the  "mass”  of  an  occupied  CA  site  is 
m,  its  energy  is  d/2mv2ltl  in  d  space  dimensions.  Then 
the  fluctuation  in  total  thermal  energy  over  a  cell  with 
n  occupied  CAs  is  v'/imi/,2,,.  (We  note  that  in  a  CA 
with  velocity  states  t  vth,  energy  fluctuations  are  pro¬ 
portional  to  density  fluctuations.)  The  corresponding 
hydrodynamic  energy  within  a  cell  of  size  tj  is  pn3v/i, 
where  vH  is  the  hydrodynamic  velocity  and  p  is  the  hy¬ 
drodynamic  density.  In  three  dimensions,  the  dissipa¬ 
tion  scale  is  tj  and  the  associated  hydrodynamic  veloci¬ 
ty  is  u,.  Also,  the  relation  between  m  and  p  is  nm 
-pr)J.  Thus,  for  thermal  fluctuations  to  be  small,  we 
must  require  that  n  »  <VRe/ A/4.  In  two  dimensions, 
the  corresponding  result  is  n  »  Nle/Al*. 

This  argument  shows  that  the  storage  and  work  re¬ 
quired  for  a  self-consistent  hydrodynamic  description 
using  CAs  is  of  order 

S-N&M*.  W  -  ,VR42/,V/7 
forCA  (2D), 

S-,VRe.  W-Ntf 
for  Navier-Siokes  (2D); 

5  -  zVRj/4/A/4,  W  - 
for  CA  (3D), 

S-.VR%4,  IF-.V’e 
for  Navier-Stokes  (3D). 

The  CA  models  approximate  fluids  that  are  by  their 
nature  necessarily  compressible.  This  means  that  an 
equation  of  state  for  pressure  is  needed.  However, 
self-consistency  requires  that  thermodynamic  pressure 
fluctuations  over  dissipation  scales  be  small.  This 
latter  condition  leads  to  results  identical  to  those  just 
obtained  by  use  of  energy  estimates.  Indeed,  it  is 
known*  that  the  rms  thermodynamic  pressure  fluctua¬ 
tions  in  a  volume  tj3  are  (pkTc/T)')1'*,  where  T  is  the 
temperature  of  the  fluid  and  c  is  the  sound  speed.  But 
the  hydrodynamic  pressure  fluctuations  over  length 
scales  of  order  rj  are  of  order  p(tv)1  ‘.  Since  */" 
-mv;n  and  c  ~  ulh,  the  previously  given  estimate-, 
apply. 

We  believe  that  these  pess.mistic  estimates  lor  high 
Reynolds  numbers  and  low  Ma^h  numbers  must  he 
overcome  belore  CAs  can  be  an  effective  modef.ne 
tool  for  complex  fluid  flows  This  can,  m  prinupic,  *'v 
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done  by  averaging  over  the  shortest  scales  <1  «  p  m 
order  to  reduce  the  number  of  degrees  of  freedom.  ’ 
However,  it  seems  that  this  renormalization  can  be 
useful  (in  the  context  of  local,  few-bit.  parallel  compu¬ 
tations)  only  if  it  does  not  generate  nonlocal,  complex 
interactions  in  the  set  of  basic  rules  defining  CAs.  Un¬ 
fortunately,  we  do  not  now  understand  why  this  kind 
of  “turbulence  transport”  modeling  should  be  either 
easier  or  more  successful  on  the  C A  lattice  than  for 
the  continuum  equations  or  for  molecular  dynamics. 

While  the  above  estimates  for  CA  simulations  of 
turbulence  are  quite  pessimistic,  there  may  be  cases  in 
which  CA  simulations  of  turbulence  may  be  effective. 
In  a  turbulent  boundary  layer,  the  local  Reynolds 
number  is  0(1)  in  the  viscous  sublayer  and  is  modest 
within  the  buffer  layer.  A  CA  model  could  be  effec¬ 
tive  in  these  regions  in  the  modeling  of  turbulent  burst 
formation  and  evolution.  However,  this  application 
requires  the  development  of  three-dimensional  CA 
models  and  suitable  techniques  to  match  the  outer  re¬ 
gions  of  the  flow. 
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I  he  recursion  relation  for  />'(')  t*  derived  readily 


rV’/dr  ♦  (v  -V)/  -  jr„V’r  (4  1) 
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I,..f„Un,c  transport  can  he  constructed  ..sin*  the  wilh  r,  =  J(l  594J  ’/C  '  *00817  I  his  relation,  which  is  usually  callc.l 

VUr‘  '*'r' ''  *■*»'  ,n  ll,,s  »*  »*T'o  hr  il.  rir.iig  an  HN.i  h..s,<l  an  algebraic  K  i  mixlel.  has  Ken  widely  used  in  lurhulcnec  modeling 

•lie  hr  IK  fur  hill,  nee  rri.xlcl  I  cl  us  assume  Dial  the  integral  se.de  /.  =  «/  f,  (I  sunder  and  Spalding.  1972;  Reynolds,  1976.  launder  el  1975)  llie 

soii.  -p.m.ls  n.  ihe  largest  llmtualing  scales  in  the  system  llns  means  that  "c.pci.mcnlally"  determined  coefficient  c,  =009  is  ijuile  close  to 

»th  >  -  SI  k  )  IS  assumed  to  he  . . .  ill  <  .1,  I  luiuii.iling  all  ..  00817  obtained  here  hy  the  RN(i  mclliml 
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S'lhslihiitnit  (*»  M  (<»  )  into  (6  2 1  gives 
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